function Analyze_pi0()
    clear all;clc;cd('./');format compact;format short g;close all hidden;

% % % %  Paras = [ r mu sigma rho gamma x0 I p  eta1 eta2 lambda T];;
    ip=11;
    dana=0.1;
    pi0=0:dana:1;
    ana=pi0;
    

    n=length(ana);

    
    Paras = Parameters;
  
 %  [xis,Xip,Phip,Ppi,Cp,Fp,Fs] 
%     Fxi0=zeros(n,1);
%     Fxi1=zeros(n,1);
    Fxis=zeros(n,1);
    Fxip=zeros(n,1);
    Fphig=zeros(n,1);
    FCs=zeros(n,1);
    FCp=zeros(n,1);
    FFp=zeros(n,1);
    FFs=zeros(n,1);
    
    Exi0=zeros(n,1);
    Exi1=zeros(n,1);
    Exis=zeros(n,1);
    Exim=zeros(n,1);
    Exip0=zeros(n,1);
    Exip1=zeros(n,1);
    Exip=zeros(n,1);
    Ephi0=zeros(n,1);
    Ephi1=zeros(n,1);
    Ephis=zeros(n,1);
    Ephip0=zeros(n,1);
    Ephip1=zeros(n,1);
    Ephig=zeros(n,1);
    ECs=zeros(n,1);
    ECp=zeros(n,1);
    EFp=zeros(n,1);
    EFs=zeros(n,1);
    
    
 
    for i=1:n
        Paras(ip)=ana(i);
         [Fxis(i),Fxip(i),Fphig(i),FCs(i),FCp(i),FFp(i),FFs(i)]=Solve_model_F(Paras);   
         [Exi0(i),Exi1(i),Exis(i),Exim(i),Exip0(i),Exip1(i),Exip(i),Ephi0(i),Ephi1(i),Ephis(i),Ephip0(i),Ephip1(i),Ephig(i),ECs(i),ECp(i),EFp(i),EFs(i)]=Solve_model_E(Paras);    
    end 
  %______________________________
  %______________________________
  %______________________________
 % [xi0,xi1,xis,xim,xip0,xip1,xip,phi0,phi1,phis,phip0,phip1,phip,Cs,Cp,pi]
 
  
 
 
  
  

          

    
     figure
    plot(pi0,Fxip,'-*b',pi0,Exip,'-sr','linewidth',3);
    xlim([0.2,1]);
    xlabel('The initial belief','Fontsize',16,'Fontname', 'Times');
    ylabel('Investment threshold','Fontsize',16,'Fontname', 'Times');
    hg = legend('PE with FGS','PE with EGS',0); optex={'fontsize', 16, 'fontname', 'Times','Interpreter','tex'};  set(hg, optex{:});
    saveas(gcf,'10-1','epsc');
    saveas(gcf,'10-1','png'); 
    
    
    
        figure
    plot(pi0,FFs,'-p',pi0,FFp,'-*b',pi0,EFs,'--k',pi0,EFp,'-sr','linewidth',3);
    xlabel('The initial belief','Fontsize',16,'Fontname', 'Times');
    ylabel('Option value','Fontsize',16,'Fontname', 'Times');
    hg = legend('SE with FGS','PE with FGS','SE with EGS','PE with EGS',0); optex={'fontsize', 16, 'fontname', 'Times','Interpreter','tex'};  set(hg, optex{:});
    saveas(gcf,'10-2','epsc');
    saveas(gcf,'10-2','png');   

    
    
end